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Abstract 

A two-dimensional atomic mass spring system is investigated for critical 
fracture loads and its crack path geometry. We rigorously prove that in the 
discrete-to-continuum limit, the minimal energy leads to a universal cleav- 
i -^h ■ age law and energy minimizers are either homogeneous elastic deformations 

or configurations that are cracked along specific crystallographic hyper- 
planes. Furthermore, we identify an effective continuum model through 
T-convergence. 
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1 Introduction 



The behavior of brittle materials is of great interest in applications as well as 
from a theoretical point of view. Such materials show an elastic response to very 
small displacements and develop cracks already at moderately large strains. In 
particular, there is typically no plastic regime in between the restorable elastic 
deformations and complete failure due to fracture. Major challenges in the ex- 
perimental sciences and theoretical studies are to identify critical loads at which 
such a body fails and to determine the geometry of crack paths that occur in the 
fractured regime. 

In variational fracture mechanics displacements and crack paths are deter- 
mined from an energy minimization principle. Following the pioneering work 
of Griffith [26], Francfort and Marigo [23] have introduced an energy functional 
comprising elastic bulk contributions for the unfractured regions of the body and 
surface terms that assign energy contributions on the crack paths comparable to 
the size of the crack of codimension one. Subsequently these models have been in- 
vestigated and extended in various directions. Among the vast body of literature 
we only mention the work of Dal Maso and Toader [20]; Francfort and Larsen 
|22j ; Dal Maso, Francfort and Toader [TH]. Determining energy minimizers of 
such functionals leads to solving a free discontinuity problem in the language of 
Ambrosio and De Giorgi [21] as the crack path, i.e., the set of discontinuity of 
the diplacement field is not pre-assigned but has to be found as a solution to the 
variational problem. In particular, these models also lead to efficient numerical 
approximation schemes, cf., e.g., pfl l6| 1281 1291 132] . 

Due to the crystalline structure of matter, under tensile boundary loads frac- 
ture typically occurs in the form of cleavage along crystallographic hyperplanes 
of the atomic lattice. On the continuum side such behavior can be modelled by 
anisotropic surface terms which are locally minimized for these crack geometries, 
see e.g. [TJ [T5J [2H]. A discrete model has been investiged by Braides, Lew and 
Ortiz [T2j, who assume that fracture can only occur in these directions and then 
calculate a limiting continuum energy: a cleavage law. This assumption leads 
to an effective one-dimenional problem which is much easier to analyze. Indeed 
in the one- dimensional setting, where discrete models describe the behavior of 
atom chains, a number of results have appeared rather recently on the literature, 
including [HI [HI DEI ED]- While by now for many atomistic models the passage 
to effective continuum models is well understood in the regime of purely elastic 
interactions, see [5j [HI [31] , not much is known on discrete-to-continuum limits 
for models allowing for fracture in more than one dimension. The farthest reach- 
ing results in that direction seem to be due to Braides and Gelli [12], who prove 
T-convergence results for scalar valued free discontinuity problems. 

However, all these ansatzes fall short of rigorous arguments that indeed in 
more than one dimension, if fracture occurs at all, then it is energetically favorable 
to cleave the specimen along particular crystallographic hyperplanes. The main 
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goal of this paper is to provide a rigorous and rather complete study of a two 
dimensional model problem. We assume that the body is a rectangular strip 
subject to uniaxial tensile boundary conditions. The atoms in their reference 
configuration shall be given by the portion of a triangular lattice in that strip 
that interact via next neighbor Lennard- Jones type potentials. This model seems 
to be the simplest model problem which (1) is frame indifferent in its vector- valued 
arguments in more than one dimension, (2) gives rise to non-degenerate elastic 
bulk terms and (3) leads to surface contributions sensitive to the crack geometry 
with competing crystallographic hyperplanes. Moreover, two-dimensional lattice 
surfaces naturally appear in the analysis of thin structures. In particular we will 
also discuss consequences of our analysis on the stability of brittle nanotubes 
under interior expansive pressure. 

Indeed we will prove that under uniaxial tension in the continuum limit the 
energy satisfies a particular cleavage law with quadratic response to small bound- 
ary displacements followed by a sharp constant cut-off beyond some critical value. 
Moreover, we will see that any sequence of minimizers converges (up to sub- 
sequences) to a homogeneous continuum deformation for subcritical boundary 
values, while it converges to a continuum deformation which is cracked along a 
crystallographic line and does not store elastic energy in the supercritical case, 
whenever the optimal crystallygraphic line is unique. The model under investi- 
gation leads, in particular, to configurations respecting the Poisson effect, which 
would not be possible in scalar models. These results justify rigorously the afore- 
mentioned assumptions in the derivation of cleavage laws as, e.g., in [13]. Finally 
we relate the sequence of discrete energy functionals to a limiting functional 
through a T-convergence result. 

The paper is organized as follows. We first introduce our discrete model and 
state our main results in Section |2j Here we already discuss different scalings of 
the boundary data and find the interesting regime where both energy contribu- 
tions, the elastic and the crack energy, are of the same order. 

In Section [3] we collect some elementary properties of the cell energy. In 
particular, we introduce a lower-bound comparison energy, called reduced energy, 
providing the optimal cell energy in dependence of the cell expansion in the space 
direction where tensile boundary conditions are imposed. 

Section [4] is devoted to the derivation of cleavage laws for the limiting min- 
imal energy. Using an elementary slicing argument we reduce the problem to 
one- dimensional segments and show that the limiting energy has a universal 
form independent of the interatomic potential. Our result is similar to the ef- 
fective one-dimesional law discussed in [13]. We obtain that the crack energy is 
anisotropic and depends explicitly on the lattice orientation. We then give finer 
estimates on the limiting energy by deriving higher order terms for the discrete 
minimal energies and show that in contrast to the limiting behavior anisotropic 
contributions also occur in the elastic regime. Moreover, our proof illustrates the 
typical behavior of brittle materials already seen in the continuum cleavage law 
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also in a discrete framework: There is essentially no plastic regime besides the 
elastic and the crack regime. More precisely, we see that for almost minimizer 
the deformation is either y/e-dose to the identity mapping (representing elastic 
response) or springs between adjacent atoms are elongated by a factor scaling 
like 4| (leading to fracture in the limit description) where e denotes the typical 
interatomic distance. In particular, here we can already see that homogeneous 
deformations or cleavage along specific lines are asymptotically optimal. 

In Section [5] we proceed to show that, under appropriate assumptions, in terms 
of suitably rescaled displacement fields indeed all discrete energy minimizers con- 
verge strongly to such continuum deformations. We provide a fine characteriza- 
tion of the crack, i.e. of the number and position of largely elongated springs. 
In the subcritical case the contribution of such springs is abitrarily small such 
that the purely elastic theory applies. For supercritical boundary values largely 
deformed springs lie in a small stripe in direction of the optimal cristallographic 
line. 

The last Section [6] finally contains our results on the limiting variational prob- 
lem. We first show that the discrete energy functionals converge to an energy 
functional defined on the continuum in the sense of T-convergence. We finally 
analyze the continuum problem under the same tensile boundary values and in 
that way we re-derive the results of Section @] and Section 

2 The model and main results 

The discrete model 

Let C denote the rotated triangular lattice 



where Rc G SO (2) is some rotation and Vi, v 2 are the lattice vectors vi = Rc&i 
and v 2 = Rc{\ G + "^ e 2), respectively. We collect the basic lattice vectors in 
the set V = {v 1; v 2 , v 2 - Vi}. The region Q = (0,/) x (0,1) C R 2 , I > 0, is 
considered the macroscopic region occupied by the body under investigation. In 
the reference configuration the positions of the specimen's atoms are given by 
the points of the scaled lattice eC that lie within Q. Here e is a small parameter 
defining the length scale of the typical interatomic distances. 

The deformations of our system are mappings y : eC n Q — > M 2 . The energy 
associated to such a deformation y is assumed to be given by nearest neighbor 
interactions as 





(1) 
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Note that the scaling factor - in the argument of W takes account of the scaling 
of the interatomic distances with e. The pair interaction potential W : [0, oo) — > 
[0, oo] is supposed to be of 'Lennard- Jones-type': 

(i) W > and W(r) = if and only if r = 1. 

(ii) W is continuous on [0, oo) and C 2 in a neighborhood of 1 with a := W"(l) > 
0. 

(iii) lim^oo W(r) = (3. 

In order to obtain fine estimates on limiting energies and configurations we will 
also consider the following stronger versions of hypotheses (ii) and (iii): 

(ii') W is continuous on [0, oo) and C 4 in a neighborhood of 1 with a := W"(l) > 
and arbitrary a' := W"'(l). 

(iii') W(r) = (3 + 0(r~ 2 ) as r -> oo, 

which is still satisfied, e.g., by the classical Lennard- Jones potential. 

In order to analyze the passage to the limit as e — > it will be useful to 
interpolate and rewrite the energy as an integral functional. Let C e be the set 
of equilateral triangles A C O of sidelength e with vertices in eC and define 
Q e = |J AgC A. By y : Q £ — > M? we denote the interpolation of y, which is affine 
on each A e C. The derivative of y is denoted by Vy, whereas we write (v)a 
for the (constant) value of the derivative on a triangle A £ C e . Then (JT]) can be 
rewritten as 

EM = ^a((v)a) + £ e boundary (y) 

(2) 

W A {Vy) dx + E £ boundary (y), 
where 

W A (F) = l - (W(|F V1 |) + W(\Fv 2 \) + W^(|F(v 2 - Vl )|)). (3) 

(Note that | A| = y / 3e 2 /4.) Here the boundary term is the sum of pair interaction 
energies \W{ ^ v ^~^ x ^ ) over nearest neighbor pairs which form the side of only 
one triangle in C e . 

Boundary values and scaling 

We are interested in the behavior of the system when applying tensile boundary 
conditions, say in ex-direction. In particular, we would like to investigate when 
and how the body breaks, i.e., 
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(1) at which value of the boundary displacement energetic minimizers are no 
longer elastic deformations but exhibit cracks and 

(2) if indeed it is most favorable for the cracks to separate the body along 
crystallographic lines. 

In order to avoid geometric artefacts, we will therefore assume that I > 4=, so 
that it is possible for the body to completely break apart along lines parallel to 
Mvi, IRv2 or R(v2 — vi) not passing through the left or right boundaries. 

Due to the discreteness of the underlying atomic lattice we have to impose the 
boundary conditions of uniaxial extension in small neighborhoods of {0} x (0, 1) 
and {1} x (0, 1), respectively: For a £ > we set 

A(a e ) = {y= (yi, y 2 ) : eC n ft -> K 2 : 

yi(x) = (1 + a e )x\ for X\ < e and X\ > I — e}. 

Note that there is some arbitrariness in this implementation of boundary values 
as one might, e.g., equally well ask that 

yi(x) = X\ for X\ < e and y\{x) — x\ + a £ l for x\ > I — e. (4) 

Such different choices will, however, not change the results of the analysis. 

Also note that there is no assumption on the second component of the bound- 
ary displacement, i.e., the atoms may 'slide along the boundary lines'. Besides 
describing a basic experiment on elastic bodies, this assumption allows for a direct 
application of our results to the stability analysis of nanotubes: 

If the rotation Rc and the length I are such that for a sequence — > the 
translated lattice SkC + (I, 0) concides with the original lattice we may view 
the system as an atomistic nanotube with macroscopic region -^S 1 x (0, 1). (Note 
that for small Ek the bending energy contributions when rolling up (0, 1) x (0, 1) 
into a cylinder are negligible as this mapping is an isometric immersion and thus 
infinitesimally rigid.) Imposing periodic boundary conditions, for arbitrary / > 
our system then models deformations of a nanotube subject to expansion of the 
diameter. 

There are two obvious choices for deformations satisfying the boundary condi- 
tions: The homogeneous elastic deformation y el (x) = {\+a £ )x and a cracked body 
deformation y cr , which, up to a boundary layer of negligible energy, is the iden- 
tity to the left and a translation by a e le\ to the right of some segment (or curve) 
passing through ft that connects a point on the lower boundary (e, I — e) x {0} 
and a point on the upper boundary (e, I — e) x {1}. It is not hard to see that 

E £ (y cl ) ~ e~ 2 W(l + a e ), E £ (y cr ) ~ e~\ 

In particular, we are interested in the most interesting regime where both of these 
energy values are of the same order, i.e., a £ is small and 

e~ 2 a 2 £ ~ e~ 2 W(l + a £ ) ~ e" 1 a £ ~ \fs. 
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In order to obtain finite and nontrivial energies in the limit e — > 0, we accordingly 
rescale E £ to £ e := eE £ . 

Conceivable alternative implementations of the boundary conditions as al- 
luded to above will then result in energy changes of order 0(e). We will account 
for all such possibilities by characterizing not only energy minimizing configura- 
tions, but more generally all configurations which are energy minimizing up to 
an error term of order O (e) . 



Limiting minimal energy and cleavage laws 

We begin our analysis with an elementary argument which yields the limiting 
minmal energy as e — > when a e j\/e — > a G [0, oo]. We first establish a lower 
bound for this energy by considering slices of the form (0, I) x {x 2 } for x 2 £ (0, 1) 
and using the reduced energy W defined by 

W(r) = mi{W A (F) : e^Fei = r}. (5) 

In a second step we show that this bound is attained. In particular, it turns 
out that the limiting minimal energy is given by elastic deformations up to some 
critical value a cr i t of the boundary displacements and by cleavage along a specific 
crystallographic line beyond a crit . 

Let 7 = max{|v 1 -e 2 |, |v 2 -e 2 |, |(v 2 — vi)-e 2 |} and v 7 G V such that 7 = |v 7 -e 2 |. 
We note that 7 takes values in [y/3/2, 1] and that v 7 is unique if 7 > y/3/2. 

Theorem 2.1 Suppose a E /\/e — > a G [0, 00]. The limiting minimal energy is 
given by 

liminf {S £ (y) : y G A(a E )} = min |-^ia 2 , ^ J . (6) 

As already motivated above, only one of the regimes is energetically favorable 
if a G {0, 00}. In the interesting case a G (0, 00) we indeed will see that in terms 
of the critical boundary displacement 



A crit 



' 2^3/3 



the limit is attained for homogeneously deformed configurations if a < a cr ; t and 
for configurations cracked along lines parallel to Rv T , if a > a cr ; t . In the special 
case that v 7 is not unique the limit is also attained if the crack takes a serrated 
course parallel to E(|, or E(-|, 

For the sake of simplicity we specialize to sequences a £ = y/ea. Without 

loss of generality we assume that R c = c ^ n /j for G [0, f ), so that 

7 = sin(</> + I ) = v 7 • e 2 . If the assumptions (ii') and (in') on W hold, we have 
the following sharp estimate on the discrete minimal energies up to error terms 
of the order of surface contributions. 
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Theorem 2.2 For e small the discrete minimal energy is given by 



[ al 9 \6a + 7a' — 2(3a — a') cos(66)]l r- % 2/3 ") . 

inf £ £ = min <^ -^a 2 + ± ^ — = ^ ^— ^-^ea 3 , — } + 0(e). 

I 27v/3 7 J 

Thus, while the zeroth order contributions in the elastic regime are isotropic, 
the higher order contributions as well as the fracture energy explicitly depend on 
the lattice orientation angle <fi. 

Detailed proofs of these results will be given in Section HI 

Limiting minimal configurations 

Our analysis of the limiting minimal energy so far showed that, depending on the 
boundary data, homogeneous deformations or completely cracked configurations 
are energy minimizing in the limit e — > 0. However, it falls short of showing that 
in fact these configurations are the only possibilities to obtain asymptotically 
optimal energies. Indeed, if v 7 is not unique, then we have already seen that 
the crack path can become geometrically much more complicated. Our next 
result shows that if v 7 is unique, energy minimizing configurations converge to 
a homogeneous continuum deformation for subcritical boundary values, while in 
the supercritical case they converge to a continuum deformation which is cracked 
along a crystallographic line and does not store elastic energy. 

The basic idea behind our reasoning will be to 'count' the number of 'broken' 
springs, i.e. the springs intersected transversally by the crack path. We see that 
the springs broken by a crack line (p, 0) + Mv 7 do not overlap in the projection 
onto the X2-axis and the length of the projection of two adjacent broken springs 
equals £7. This leads to a fracture energy of approximately — . If we assume that 
the cleavage is not parallel to Mv 7 we conclude that some springs in v 7 direction 
must be broken, too. If we consider the adjacent triangles of such a spring and 
their neighbors we find that the projection onto the avaxis of broken springs 
overlap. A careful analysis of this phenomenon then shows that every broken 
spring in v 7 direction 'costs' an additional energy of « 2ef3^^-, where P("f) is 
the geometrical factor 

P(7) = V3^) . (7) 

(Note that P (7) = if and only if 7 = ^ in accordance to the above considera- 
tions.) 

In order to give a precise meaning to the convergence of discrete to continuum 
deformations, to each discrete deformation y : sC M 2 we assign - as mentioned 
above - the affine interpolation y on each triangle A G C e . Accordingly, to the 
rescaled discrete displacements u : eC — > M 2 with y = id + ^/eu (id denoting the 
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identity mapping id(x) = x) we define u to be its affine interpolation on each 
triangle A G C £ . 

In the cracked regime we may of course only hope for a unique limiting de- 
formation up to translation of the crack path. However, without an additional 
mild extra assumption on the admissible discrete configurations or their energy 
even this cannot hold true, as apart from the crack, parts of the specimen could 
flip their orientation and fold onto other parts on the body at zero energy. In 
order to avoid such unphysical behavior we add a frame indifferent penalty term 
X > to Wa with x > c x in a neighborhood of 0(2) \ 5*0(2) and x = in a 
neighborhood of 5*0(2): 

W AtX = W A + x- (8) 

We set 

8f(y) = -L- [ W A , x (Vy) dx + eE^ d ^(y), 

for u G A(a £ ). More generally than a sequence of minimizers we will consider 
sequences (y e ) of almost minimizers that satisfy 

£?(y.) = inf{£?(y) : y G A(a £ )} + 0(e). (9) 
For those deformations we will show in Section [5j 

Theorem 2.3 Assume that W satisfies (i), and (Hi'). Let v 7 be unique, 
a e = yfea, a ^ a crit and suppose (y e ) satisfies fl9]). Let u e such that y e = id+y^i^. 
Then there exist u £ : Q — > 1R 2 with \{x G Q £ : u E (x) ^ u e (x)}\ = 0(e) such that: 

(i) If a < a cr it, then there is a sequence s e G K. such that 

\\u £ -(0,s £ )-F a -\\ m[n) ^0, 

'a ' 



where F a 



-f 



(ii) If a > a cr i t , i/ien t/iere exzst sequences p £ G (0,/), s e ,t e G 1R snc/i i/iai 
(Ps)0) + Mv 7 intersects both the segments (0,1) x {0} and (0,1) x {1} and, 
for the parts to the left and right of (p £ , 0) + Rv 7 

:= {x G Q : < X\ < p £ + (v 7 • e\)x2} and 

:= {x G : p E + (v 7 ■ ei)x2 < Xi < /} , 

respectively, we have - possibly after rotating y £ by tt on or - 

\\u £ - (0,^)11^1(0(1)) + \\u £ - (al,t £ )\\ H i {ni 2 )] 0. 

Note that a rotation by ix on nw, j = i,2, might be necessary as on each of 
these sets there are, up to translation in ^-direction and a small correction in 
a boundary layer of negligible energy, two deformations respecting the boundary 
conditions that do not store elastic energy: y(x) = x and y(x) = —x in 
respectively, y(x) = x + a E le\ and y(x) = —x + (2 + a £ )le\ in 
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The limiting variational problem 

We finally address the more general question if not only the minimal values 
or the minimizers but the whole energy functionals converge to a continuum 
energy functional in a variational sense. Furthermore, we analyze the limiting 
problem independent of its discrete approximations. The results announced here 
are proved in Section [6j 

Our convergence analysis applies to discrete deformations which may elongate 
a number scaling with ^ of springs very largely, leading to cracks of finite length in 
the continuum limit. On triangles not adjacent to such essentially broken springs, 
the defomations are y^-close to the identity mapping, so that the accordingly 
rescaled displacements are of bounded L 2 -norm. Note that the first of these 
assumptions can be inferred from suitable energy bounds. By way of example, 
however, we will see that this cannot be true for the displacement estimates in 
the bulk: The sequence of functionals (S £ ) is not equicoercive. Nevertheless, 
it is interesting to investigate this regime in order to identify a corresponding 
continuum functional which describes the system in the realm of Griffith models 
with linearized elasticity. As a particular case of Theorem 12.31 we mention that 
(£ e ) is still mildly equicoercive. 

Recall that the space SBV(Q; IR 2 ), abbreviated as SBV(Q) hereafter, of spe- 
cial functions of bounded variation consists of functions u G L 1 (fi;IR 2 ) whose 
distributional derivative Du is a finite Radon measure, which splits into an ab- 
solutely continuous part with density Vw with respect to Lebesgue measure and 
a singular part D^u whose Cantor part vanishes and thus is of the form 

D j u = [u] ® v u U l [J u , 

where H 1 denotes the one- dimensional Hausdorff measure, J u (the 'crack path') 
is an T^-rectifiable set in Q, v u is a normal of J u and [u] = u + — u~ (the 'crack 
opening') with u being the one-sided limits of u at J u . If in addition Vw G L 2 (Q) 
and l-L l {J u ) < oo, we write u G SBV 2 (VL). See [3] for the basic properties of these 
function spaces. 

The sense in which discrete displacements are considered convergent to a 
limiting displacement in SBV is made precise in the following definition. 

Definition 2.4 Suppose u £ is a sequence of discrete displacements such that the 
corresponding deformations y £ = id + yfeu £ are uniformly bounded in L°° . We 
say that u £ converges to some u G SBV 2 (Q): u £ — > u, if 

(i) Xn E u £ in L 1 ^) 

and there exists a sequence C* C C £ with #C* < j for a constant C independent 
of e such that 

(H) ||Vu £ |U2 ( n e \u AeC| A) < C. 
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The main idea will be to separate the energy into elastic and crack surface 
contributions by introducing a threshold such that triangles A with (y)& beyond 
that threshold are considered as cracked and y is modified there to a discontinuous 
function. 

Consider the limiting functional 

S(u) = -j=J^Q(e(u))dx + j j Y^^-^dH 1 (10) 

for u G SBV 2 (Q), where e{u) = | (Vm t + V«) denotes the symmetric part of 
the gradient. Q is the linearization of about the identity matrix Id (see 
Section H] for its explicit form). For the sake of simplicity we again suppose that 
a £ = -y/ea for all e. We then have the following T-convergence result: 

Theorem 2.5 Let a £ = \fea and a G [0, oo). 

(%) If (u £ ) is a sequence of discrete displacements such that y £ = id + \feu £ G 
A(a £ ) and u £ — >• u G SBV 2 (Q) with ui(0, •) = and «i(/,-) = la (in the 
sense of traces), then 

liminf£ e (u e ) > £{u). 

(ii) For every u G SBV 2 (Q) with ui(0, ■) = and ui(l,-) = la (in the trace 
sense) there is a sequence {u £ ) of discrete displacements such that y £ = 
id + ^Jeu £ G A(a £ ), u £ u G SBV 2 (Q) and 

\xm£ £ (u £ ) = S{u). 



We note that the sequence (S £ ) is mildly equicoercive in the sense that low energy 
sequences (satisfying admit a subsequence converging in the sense of Defi- 
nition [231 by Theorem 12.31 (the convergence is even stronger in this case). Due 
to the frame indifference of W, (S £ ) is not equicoercive as the following example 
shows. 

Example. Assume that the specimen satisfying the boundary conditions is bro- 
ken into three parts by two even cracks where the middle part is subject to a 
rotation R ^ Id so that 



Vy £ = R for p < X\ < q, < p < q < I. 



In particular, the energy of the configuration is of order 1. But for p < X\ < q 

1 



WuJx) 



(R - Id) 



— > oo for e — > 0. 



Thus, Vu £ is not bounded in L 1 and so u £ does not converge. 
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Finally, the limiting functional 8 can also be analyzed directly without re- 
course to the approximating functionals. We determine the minimizers and prove 
uniqueness up to translation of the specimen and the crack line for the boundary 
conditions 

Ui = for X\ — and U\ = al for X\ = I. (11) 

Theorem 2.6 Let v 7 be unique and a ^ a crit . Then min£ = min j^=a 2 , 
and all minimizers of 8 subject to ( TTTT) are of the following form: 

(i) If a < a crit , then 

for some s G K. 

(ii) If a > a cr it, then 




(0, s) for x to the left of (p, 0) + Rv 7 , 
(aZ, t) for x to the right of (p, 0) + Rv. 



/or some s,i G R and p G (0,/) such that (p, 0) + Mv 7 intersects both the 
segments (0,1) x {0} and (0,/) x {1}. 

An analogous result for nonlinear but isotropic energy functionals has been ob- 
tained recently by Mora-Corral [27]. 

We close this introductory chapter emphasizing that all the optimal configura- 
tions found in Theorem 12 . 3 1 and Theorem 12. 6 1 by minimizing the energy without a 
priori assumptions show purely elastic behavior in the subcritical case and com- 
plete fracture in the supercritical regime. In particular, the elastic minimizer 
u cl shows elongation a in ei-direction and compression — | in the perpendic- 
ular e 2 -direction, a manifestation of the Poisson effect (with Poisson ratio |), 
which cannot be derived in scalar valued models. On the other hand, the crack 
minimizer u cr is broken parallel to Mv 7 which proves that cleavage occurs along 
crystallographic lines. 

3 Elementary properties of the cell energy 

We collect some properties of the cell energy Wa and the reduced energy defined 
in for W satisfying the assumptions (i), (ii) and (iii). 

Lemma 3.1 Wa is 

(i) frame indifferent: W A (QF) = W A (F) for all F G M 2x2 ; Q G SO{2), 



u e \x) = (0, s) + F a x 
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(ii) non-negative and satisfies W A (F) = if and only if F £ 0(2) and 
(Hi) liminfi^i^oo W A {F) = (5. 

Proof, (i) is clear. For (ii) it suffices to note that v T F T Fv = 1 for three vectors 
v, no two of which are collinear, implies that F T F = Id. (iii) can be seen by 
noting that if \F\ — > oo, then for at least two vectors v £ V one has \Fv\ — > oo. 

□ 

We compute the linearization about the identity matrix Id: 
Lemma 3.2 Let F = Id + G for G £ R 2x2 . Then for \G\ small 

W A (F)= 1 -Q(G) + o(\G\ 2 ), 

where Q(G) = ff (3g 2 n + 3g 2 2 + 2g u g 22 + 4 (m^if 

In particular, Q{G) only depends on the symmetric part (G T + G) /2 ofG. Q 
is positive semidefinite and thus convex on M 2x2 and positive definite and strictly 
convex on the subspace M 2 *^ °f symmetric matrices. 

Proof. Let v £ V and G £ IR 2x2 small. We Taylor expand the contributions 
W(|-Fv|) to the energy W A : 



W(\{I& + G)v\) = W V(v, (Id + G T )(Id + G)v) 



W"(l) I G T + G 



2 

Now using the elementary identity 



2 



o(\G\ 



(v l7 H^) 2 + (v 2 , Hw 2 ) 2 + ((v 2 - Vl ), H(v 2 - Vl )) 2 

3 (12) 
= - (2trace(i/ 2 ) + (traced) 2 ) 
8 

for any symmetric matrix H £ M 2x2 , we obtain by summing over v £ V 
W A (F) = _._._.[ 2trace ( ( J +( trace ) ) + o(\G\ 



2 2 
1 -Q(G) + o(\G\ 2 ). 



As Q(G) > ff (2g 2 ± + 2g 2 2 + (g l2 + g 2 i) 2 ), Q is positive semidefinite on M 2x2 and 
positive definite on ffi 2 ^- D 

As a consequence, we have the following properties of the reduced energy W. 
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Lemma 3.3 The reduced energy satisfies 
(%) W(r) = |r| < 1. 

(ii) For r > 1 one has 

W(r) = W A A)) + o((r - I) 2 ) = J(r - l) 2 + o((r - l) 2 ). 

fmj limir^oo W(r) = /?. 

Proof, (i) If r < 1, then one can choose Q G 50(2) with ejQei = r and 
so < I¥(r) < W A (Q) = 0. If \r\ > 1, then W(r) > for otherwise there 
would be a sequence G M 2x2 with efF^ei = r and Wa(-Pa:) — ► 0. But then 
dist(Ffc, 0(2)) by (ii) and (hi) of Lemma [37TI and thus, up to subsequences, 
Fk — > F G 0(2) with ejFe\ = r, which is impossible. 

(ii) This discussion shows that in fact for any 6 > there exists n > such 
that W/\(F) > 5 whenever dist(F, 0(2)) > rj. Now since W(r) — > as r \ 1, we 
obtain that, for sufficiently small r > 1 and 5 > 0, any F with W/^-F 1 ) < VF(r) + 5 
is contained in a small neighborhood of 0(2). If in addition efFei = r holds, then 

in fact, F must be close to Id or to P = . In particular, by continuity 

of W, the infimum on the right hand side in the definition of W is attained for 
those r. 

We now fix such an r > 1 near 1 and choose F = Id + G such that W(r) = 
W A (F) and ejFei = r. As W A is invariant under the reflection P, we may 
without loss of generality assume that G is small. Then Lemma 13.21 yields 

W A (F) = | + 3, 2 2 + 2g u g 22 + 4 {^^-)^ + o(|G| 2 ). 

We find that gn — r — 1, g i2 + g 2 i = °( r — 1) an d #22 = — ^dn + o(r — 1) and F 
satisfies 

-F T + F (r 



V 



o(r — 1, 



W A (F) = W A [ FT ^ F ) + o((r - l) 2 ) 



with energy 



= ^(r-l) 2 + o((r-l) 2 ). 

(hi) This is immediate from Lemma [3.1 ( hi). □ 
Under strengthened hypotheses on W we have the following expansion: 
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Lemma 3.4 IfW in addition satisfies the assumptions (ii 7 ) and (Hi'), then for 
r > 1 close to 1 we have 

W(r) = a(r ~ 1)2 + -L (6a + 7a' - 2(3a - a') cos(60)) (r - l) 3 + 0((r - l) 4 ), 
4 108 V / 



cos cp — sin t 
sin <b cos (A 



where <f> is such that ~- 

Proof. Let s = r — 1. By definition, 

W(r) = mm {W A (F(s,x,y,z)) : x, y, z G M.} , 

where F(s,x,y, z) =( 1 ^ s f!u^ )■ ^ ue ^° the quadratic energy growth near 

S0(2), we need to minimize only over x,y,z with \z\, \/s\y\ < Cs for a 
constant C large enough. Indeed, as W&(F(s, 0, 0, 0)) = 0(s 2 ), for a minimizer 
one has without loss of generality dist(-F(s, x, y, z), SO (2)) = O(s). But then 
a/(1 + s) 2 + (z ± y) 2 = 1 + O(s), which implies |z ±y\ = 0(y/s) and so \z\, \y\ = 
0(y/s), and also \/ (I + x) 2 + (z ± y) 2 = l + 0(s), which then implies ±(l + x) = 
1 + 0(s) and thus without loss of generality x = O(s). Finally using that the 
scalar product (1 + s)(z + y) + (1 + x)(z — y) = 2z + 0(s 3 ^ 2 ) of the two columns of 
F(s, x, y, z) in absolute value is also bounded by 0(s), we obtain that \z\ = O(s). 

Set x = — | + sxi, y = y/syi, z = sz± with \xi\, \y±\, \z\\ < C. Explicit 
calculation gives 

W A (F(s, x, y, z)) = ^ (8 + 3x1 + %vl + 122? + 6(xi + ylf) s 2 + 0(s 3 ). 

Since a > 0, we thus obtain that this expression is minimized in x±,yi,zi with 
x 2 ,yl,zf = O(s) and we may set x\ = \fsx2, y\ = \/sy2 and z\ = \fsz2 with 
1^2 1 > 1 2/2 1, 1 -2^2 1 < C for some C > 0. Explicit expansion in powers of s then yields 

W A (F(s,a;,y,z)) 
as 2 / 

= + (48a + 56a' - 16(3a - a') cos(60) 

4 864 V v 1 v y) 

+ 3a (81x| + 72y 2 + 108z|) )s 3 
+ — ( (9ay 2 + a' + (3a — a') cos(60)) x 2 



+ 2(3a - a') sin(60)z 2 ) ) s 7/2 + 0{s^ 



as 2 



+ 77^-f6a + 7a' — 2(3a — a') 008(60)") s 3 
108 V / 



108 

9° / 2 \ 3 a l/2 s3 3a / 



+ — (xl + 2Av^x 2 J s 3 + + — [4 + 2B^z 2 ) s 3 + 0(s 
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for A and B bounded uniformly in s and so 

W A {F(s,x,y, z)) 

= ~ + 108 V a + 7a ' ~ 2<y3a ~ ^ cos ( Q ^J s 
+ |(x 2 + A*)** + ^ + f (* + B^) V + 0(A 

Minimizing with respect to 0:2,2/2 and ^2 we finally obtain that 

W(l + s) = ^- + -L-(ea + 7a' -2(3a-a')cos(60))s 3 + O(s 4 ). 
4 108 V / 

□ 

The following lemma provides useful lower bounds for the energy Wa and the 
reduced energy W. 

Lemma 3.5 For allT > 1 one has: 

(i) There exists some c > such that cdist 2 (F, 0(2)) < W A {F) for all F e 
R 2x2 satisfying \F\ < T. 

(ii) For 5 > small enough, there is a convex function V > with V(r) < W(r) 
forr < T and such that the second derivative V+(l) from the right at 1 exists 
and satisfies V£(l) = f - 26. 

(Hi) If in addition W satisfies assumptions (ii') and (Hi'), then there exists a 
convex function V > with V(r) < W(r) < V(r) + 0((r - l) 4 ) for r <T. 

(iv) For p > there is an increasing, subadditive function ip p : [0, 00) — > (0, 00) 
which satisfies ip p {r) — p < W(r + 1) for all r > and ip{r) = (3 for all 
r > c p for some constant c p only depending on p. 

Proof, (i) Let F G M 2x2 satisfying \F\ < T. By polar decomposition we find 
R G 0(2) and U = V F T F symmetric and positiv definite such that F = RU. A 
short computation yields \U — Id| = dist(F, 0(2)). Assume first \U — Id| < rj for 
77 > small enough. Since Wa{F) is invariant under rotation and reflection we 
obtain applying Lemma 13.21 

W A (F) = Wa{R t RU) > ^Q(U - Id) + o(\U - Id| 2 ). 

Noting that Q grows quadratically on M 2 *^ (see Lemma |3T2|) we obtain a constant 
c\ > such that for \U — Id| < 77 

Wa(F) >cx\U- Id| 2 = ci dist 2 (F, 0(2)). 
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Consider the compact set M := {F G M 2x2 , dist(F, 0(2)) > rj, \F\ < T}. W A 
attains its minimum on M, which is strictly positiv by Lemma I3.1( ii). This 
provides a second constant C2 > such that for all F G M 

W A (F) >c 2 \U- Id| 2 = c 2 dist 2 (F, 0(2)). 

Taking c = min{ci, C2} yields the claim. 

(ii) We construct such a function directly applying Lemma 13.31 

{0 for r < 1, 

(f-5)(r-l) 2 forl<r<l + 7/, 

(f - 5) r] (2r - 2 - rj) for r > 1 + 77, 

when 77 > is sufficiently small. 

(iii) With/(r) := ^^^ + T ^(6«+7« / -2(3a-a / )cos(60))(r-l) 3 -O(r-l) 4 
for sufficiently large O, Lemma 13.41 shows that we can choose 

for r < 1, 

V(r) = { f(r) for 1 < r < 1 + rj, 

f{\ + 7]) + 77) (r - 1 - r/) forr>l + ?7, 

when 77 > is sufficiently small. 

(iv) We define 

f]r for < r < 



for r > 



'/ 



for some rj > (depending on p) such that ^ — p < W. Then we set ip p (r) = 
ip(r + 1). As ip p is a concave function with ip p (0) > 0, it is subadditive. □ 



4 Limiting minimal energy and cleavage laws 

We now prove Theorems 12.11 and 12.21 on cleavage laws and fine energy estimates. 
Limiting minimal energy 

We can classify (or 'color') all triangles in C e into two types, say 'type one' and 
'type two', such that all triangles of the same type are translates of each other. 
Then only triangles of different type can share a common side. Denote the sets 
by Ce^ and , respectively. 

Proof of Theorem \2.1\ We first show that the expression on the right hand side is 
a lower bound for the limiting minimal energy. For every deformation y G A(a £ ) 
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we have by (j2D and (J3]) 



g £ (y)>^L [ W A (Vy)dx. 

V 3£ ^O e n(0,0x(e,l-e) 

Let < 6 < j and choose R so large that W(r) > (3 - 6 if r > R. Define 
to be the set of those triangles A of type one for which at least one side in the 
deformed configuration y(A) is larger than 2Re. By I C (e, 1 — e) we denote the 
set of those points x 2 for which there exists x\ G (0,/) such that (xi,x 2 ) lies in 
one of these triangles. 

We can then estimate the energy integral by splitting the ^-integration into 
a first part where x 2 ^ I and a second part with x 2 G I. 

1. If x 2 I, then all sidelengths of y(A) for a triangle A whose interior 
intersects the segment (0, /) x {x 2 } are less or equal to 4Re. This is clear for 
triangles of type one by construction. For triangles of type two it follows from 
the fact that the two sides of A intersecting (0, /) x {£2} ar e also sides of triangles 
of type one and therefore bounded by 2Re. The third side is thus less than 4Re, 
too. 

It is elementary to see that for F G M 2x2 

\e^Fei\<8R, if |v T Fv| < 4R for all v G V. (13) 

Indeed, if Ai, A2 are the eigenvalues of |(F T +F), then by (TT2"]) one has f (A^+Ai 2 .) = 

I trace (±(F T + F)) 2 < 3 • (4R) 2 and thus |efFei| < max{|Ai|, |A 2 |} < 8R. 
Consequently, for almost every x 2 £ I we have eJ'Vy(xi, x 2 )ei < 8R for all 
Xi6(0,Z). 

By Lemma l3.5( ii) choose a convex function with V(r) < W(r) for r < 8R 
and V^'(l) = f -25. For x 2 G (e, 1-e) define ttf C (0, Z) such that fif 2 x {x 2 } = 
Q e fl (0, /) x {22}. Then for the first part one obtains, if a < 00, by convexity of 
V 

-4- / / W A Ny)dx l dx 2 >^ f [ V(e T l Vye 1 )dx 1 dx 2 

V3e J(e,i-e)\i JnP V3e J(e,i-e)\i Jn x e 2 



>-i- / |fif|y(l + a £ )& 2 

2 



V3e 
2 

7! 



(e,l-e)\/ 

(l-2,-|/|)(Z-25)(n'(l)a 2 + (e)) 



-> -^=(1 - \i\)iv"(iy 



(14) 



as £ — > 0. It is not hard to see that this asymptotic estimate remains true also 
for a = 00. 
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2. On the other hand, the energy of the second part can be estimated by the 
energy of all springs lying on the side of a triangle in Ce , which yields 



[ [ W A (Vy)dx 1 dx 2 >2(/3-5)e#CP, 



(15) 



as the length of at least two springs in each of these triangles is larger than Re in 
the deformed configuration. Now the projection of any triangle onto the avaxis 



is an interval of length ej, and so £7#Q > |/|, i.e. 
4 



)iection ol 
(l1 > |/|, i- 

W A (Vy) dx 1 dx 2 >2((3 -d)^ 1 ]!]. (16) 



V3e J i Jnp 
Summarizing ( Fl4l) and ( Fl6l) we find 

liminf mm{£ e (y) : y G A(a e )} 

E—tOO 

> min {4= (f " 25) la\l - |/|) + 209 - ^"Vl : l J l G [0, 1] 



Now 5 — )■ shows 



al j 2/3 



liminf min{£ £ (?/) : y G „4(a e )} > min < —=a , — > . 

lV3 7 J 

This establishes the lower bound. 

It remains to prove that the right hand side in Theorem [2J] is attained for some 
sequence of deformations. In order to do so, we consider two specific sequences 
of deformations. First, for a < oo let 



yf(x) = (ld + F^)x = f 1 * " 
By Lemma [3T3lf ii) we have that W A (F) = ja 2 £ + o(e) and so 



(17) 



Yim£ £ (yf) = ^=a 2 

by ©. 

To define y cr we choose any line (s, 0) + IRv 7 intersecting both the segments 
(0,1) x {0} and (0,/) x {1} (as in Theorem 12.6ft . This is possible since I > -j=. 
Let a > and set 



. x for ce to the left of (s, 0) + Rv~, , 

yf(x) = \ (U 
x + ajei for x to the right of (s, 0) + IRv 7 
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for atoms x with e < x\ < I — e. Except for a negligible contribution from the 
boundary layers, the energy of this configuration can be estimated as in Step 2 
of the proof of the lower bound: It is given by the energy of springs intersecting 
(s, 0) + Rv 7 , i.e., by the two springs lying on the boundary of the triangles of 
type one which are intersected by (s, 0) + Rv 7 . These springs are elongated by a 
factor scaling with a e /e, thus yielding a contribution /? in the limit e — > 0. □ 

Fine estimates on the limiting minimal energy 

Assume now that W in addition satisfies assumptions (ii') and (iii'). In order to 
investigate a deformation y again we let C £ and denote the set of triangles 
A (of type one respectively) for which at least one side in y(A) is larger than 
2Re, where now the threshold value R > 1 is chosen in such a way that cr := 
infjW(r) : r > R] > |. According to Lemma [3.5( iii) we may choose a convex 
function V such that 

< V(r) < W(r) < V{r) + 0((r - l) 4 ) for r < SR. (19) 

As in (|T3|) we observe that |ef (y)A&i\ is bounded by 8R on triangles with bond 
length not exceeding 4Re and thus lies in the convex regime of V. Moreover, we 
find that every triangle in C £ provides at least the energy f A W&(Vy) > crE. 

For given < 77 < a we also define R £>r) = as a threshold for triangles we 
consider 'essentially broken': 

C 6 ,„ = {A G C e , |Vy e v| > R eyV for at least two v G V} . (20) 

The minimal energy contribution of all the springs on such a triangle in C £j7? is 
given by 

2/3% = 2inf |V(r) : r > ^7^} e = (2/3 + 0{e))e 

by the assumption (iii') on W. By / C (e, 1 — e) we denote the set of points x 2 
for which the segment (0, /) x {^2} intersects a broken triangle (of type one) in 
Ce ■ In addition, we say %2 G J* 7 C / if one of the intersected triangles lies in 

With these preparations we can now proceed to prove Theorem 12.21 

Proof of Theorem \2.2l Let S £ (y) = inf £ e + 0(e). Inspired by ( !T4|) and (|T5|) we 
establish a lower bound for the energies additionally taking the set I \ I v into 
account. Since the sidelength of any triangle whose interior intersects (0, 1) x (J\ 
J 77 ) is bounded by 4_R £j?? , we find 

|efVy(xi,x 2 ) e x | < 8R ejTj 



20 



for all (xi,X2) G (0,/) x ( I \ I v ) as in (TT3T) . Let k = k(x 2 ) count the number of 
triangles in C e on the slice (0,/) x {x 2 }, x 2 G I \ I v and define C* 2 C (0,/) such 
that ((0,/) x {x 2 }) D LheCe A = x M- Then 

/ efV^f(xi,a; 2 )ei > (1 + a/so) (Z + 0(e)) - SkeR £tV 
Jn* 2 \cT 2 



and so 



l + Vifa- 8 *^ ^ + Q(Vi) ) ) I 



A convexity argument as in the proof of Theorem 12. II on slices (0, /) x {x 2 } with 
x 2 E (e,l — e) \ I and on the unbroken part {VL* 2 \ C* 2 ) x {x 2 } of slices with x 2 
in / \ I v then shows that 

S E (y) > A{l ~ 2£) V(l + y/ia){l -2e-\I\) + G^I \ P\ + + 0(e), 

y/3e 7 

(21) 

where 

= mm -^V 1 + U S — - + 0(y/e) + 



,y/3e V V 1 J J 7 

We note that this minimum exists and can be taken over 1 < k < Kq for some 
Kq g N large enough and independent of rj as — — > oo for fc — > oo. We choose 
< 77 < a large enough such that 



/a 2 I al f 8k. \ 2 /ccr 
— =a < min —= a —la — n) \ H 

v/3 i<fc<^o V i J 7 

Recalling that, by flT9]) and LemmaES! -^/(l+^/er) = ^W{l+y/er)+0(e) -> 
-^r 2 uniformly in r on bounded sets in M., we see that thus G Vt£ exceeds the elastic 
term -^LV(1 + y/ea) for e sufficiently small. So from ( 121]) we obtain 

4/ 2/^ 
£(!/) > -j=-V(l + >/ia)(l - 2e - |/"|) + -?-\r>\ + O(e). (22) 
V3e 7 

As -^^(1 + s/ea) -> -^=a 2 and f3 v ^ f3 for all 77 > 0, for e small enough we 

thus obtain inf£ £ > *Lv(l + y/Ea){l-2e) + 0(e) = 4^\W(l + ^ea) + 0(e) or 

inf £ e > 2/3+ ^ £ - > (1 — 2e) = — + O(e), respectively, depending on a. 

Applying flEJ) and (PSj) we then get indeed inf £ £ = + y/ea) + 0(e) 

or inf £ e = — + O(e), respectively. The claim now follows from Lemma [3.41 . □ 

Remark. From the proof of Theorem 12.21 especially taking ( 1P7|) and (I18p into 
account, it follows that Theorem 12.21 still holds if £ £ is replaced by 
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5 Limiting minimal energy configurations 

Throughout this section we will assume that a £ = y/ea, y £ is a sequence of 
deformations satisfying (jUJ), the threshold value R is chosen as above Equation 
( |T9|) and that C £ is defined accordingly. 

For a rescaled displacement u we denote by C (0, 1) for \x > the 
set of x 2 such that there is precisely one triangle A X2 e Ce with int(A X2 ) D 
((0,/) x {x 2 }) ^ and 

/ e x Vu(x\, x 2 )ei dx\ < Iji. (23) 
Note that .D^ 1 C I v for small enough: For x 2 € -D M we have 

/ ef Vy(xi, x 2 )e 1 dx\ > y/el(a — fi) + 0(e) 
Jc* 2 

and using the arguments in ( fl3l) we see that for given rj we can choose small 
enough such that A X2 e C efl and thus x 2 E I v . We also define C^„ C C £)7? as the 
set of those essentially broken triangles A for which there exists some x 2 G 
such that int (A) fl ((0,/) x {£2}) 7^ 0- The projection of a triangle A onto the 
linear subspace spanned by v^r is an interval of length ^e. We denote the center 
of this interval by m a • 

The following lemma gives sharp estimates on the number of broken triangles 
and their position. 

Lemma 5.1 Let v 7 be unique and a 7^ ct C rit- Let u £ be a minimizing sequence 
satisfying 

E £ (\d + y/eu £ ) = inf £ £ + 0(e). 
(i) If a < a CTit , then e#C £ = 0(e). 

(ii) If a > a cr i t; then \I V \ = 1 — 0(e) for < 77 < a. Furthermore, for fi 
sufficiently small, ej^ (C £ \ C£„) = 0(e) and 

sup{|m Al -m A2 |,Ai,A 2 e C^ v } = 0(e). 

Proof, (i) Using ( l2T|) we find 
4/ ~ 

S £ (y £ ) = -j=-W(l + s/la) + 0(e) 

> A{l ~ 2£) W(l + V£a)(l -2e - |/|) + min (g„ £ , —} \I\ + 0(e) 
V3e I 7 J 

= + y/ia)(l ~ + min <^ G ve , — } \I\ + 0(e). 

V3e I 7 J 
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An elementary computation yields, whenever e is small enough, 



|/| < (min/c^,^! - -^-W(l + Via)) ■ 0{e) 
= Lin |g„, £ , ^} - ^=a 2 + o(l)) 1 ■ 0(e) = 0(e), 

(The argument leading to fl22|) together with a < a crit shows that the term in 
parentheses is bounded from below by a positive constant independent of e). 
Then the elastic energy is -j^W(\ + \^ea) + 0(e) and consequently, the crack 
energy coming from triangles in C £ is of order 0(e). As every broken triangle in 
C £ provides at least energy ecu we conclude £#C| = 0(e). But then, possibly 

—(2) 

after replacing R by 2R, also ej^C £ = 0(e) as those triangles are neighbors of 
broken triangles of type 1. 

(ii) Using (122j) we find after without loss of generality choosing 77 sufficiently 
large 

913 4/ ~ 2fi r > 

£e(Ve) = j + 0(e) > ^^(1 + v^)(l - 2e - \P\) + F\. 

So for e small enough we obtain 

al 2 , 2/3V 1 



1 - \P\ < l-j=a 2 + o(l) - ^ j ■ 0(e) = 0(e) 

since a > a cri t- Consequently, the crack energy from triangles in C £)V is given by 
— + 0(e) and thus the energy contribution from C £ \ C £tV is of order 0(e). As in 
(i) we find e# (C £ \C £ , V ) = 0(e). A convexity argument yields that the energy 
of a slice in I v \ is larger or equal to 



— + mm\c R ,^W(l + y/epL) \ + 0(e) = — +mm\c R , ^=^ 2 \ +o(l). 
7 I V3e J 7 I V3 J 



It follows |7' 7 \Z) M | = 0(e). We conclude |£)^| = 1 — 0(e), whence the crack energy 
from triangles in C£ v is given by — + 0(e) and then also e# (C £ \ C£ v ) = 0(e). 

Finally, we concern ourselves with the projected distance of triangles in . 
We first note that it suffices to show 

sup{|m Al -m Aa | : A l5 A 2 G n C« } = 0(e) 

since for a suitable 77 > rj for any A 6 D Ci 2 '' there is a A G H Ci 1 ^ with 
|m A — ?ri A I < £■ Let x 2 , z 2 G ^2 < ^2 with z 2 — ^2 < Oe and |m Al — I > 
for the corresponding broken triangles Ai, A 2 G C^. We may assume if a triangle 
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intersects (0,/) x {z 2 } or (0,/) x {x 2 } then its interior does so, too. Denote by 
d = 7~ 1 |?tia 1 —rri/x 2 | the distances of the centers in v 7 -projection onto the Xi-axis. 

Let x\,z\ G (0,/) be the points on the slices (0,/) x {x 2 } and (0,1) x {z 2 } 
satisfying ir v ±(xi,x 2 ) = and ir v ±(zi,Z2) = mA 2 , respectively, where 7r v ± 

denotes the orthogonal projection onto the linear subspace spanned by vrr . Let 
w — ei ■ v 7 |x 2 — z 2 \/j- Then the v 7 -projection of z = (z\, z 2 ) onto the x-slice is 
given by (zi, x 2 ) with z\ — Z\ — w. Then d — \%i_ — Z\\ and without restriction we 
may assume X\ > Z\. 

Let s e = We now consider the area bounded by the parallelogram with 

corners (z\ + s E ,x 2 ), (x\ — s e ,x 2 ), (z± + d — s e ,z 2 ), (zi + s e ,z 2 ). It is covered 
by — 1 stripes of width in v 7 -direction consisting of lattice triangles 
intersecting the parallelogram, the first of these stripes touching A 1; the last 
one touching A 2 (note that if 7J = the parallelogram is degenerated to a 
segment). For the intermediate stripes (1251) shows that 

Vi (t, x 2 ) <t+ y/elfi Wt < x\ — s £ and 
Vi(t, z 2 ) >t + y/el(a — fi) Wt > zi + s £ . 

This shows that if (t, x 2 ) and (t + w, z 2 ), X\ — d + s e < t < x x — s £ lie in the bottom 
and top triangles of some intermediate stripe, respectively, which are unbroken 
by construction of D^, then 

\y(t + w,z 2 ) -y(t,x 2 )\ > yi{t + w,z 2 ) -yi(t,x 2 ) > w + y/el{a - 2fi) ~ ^fe. 

Consider the atomic chains in v 7 direction that lie on the boundary of these 
stripes. They are of length r )~ 1 {z 2 — x 2 ) + 0(e) < Ce y/e. So there is a constant 
c > such that each of these chains contains at least one spring elongated by 
a factor of more than By passing, if necessary, to a lower threshold fj > rj, 
we obtain that the triangles sharing such a spring are broken and additionally 
one neighbor of each. As broken triangles for such springs on neighboring chains 
might overlap, we only consider every second atom chain and denote the set of 
type one triangles adjacent to such a spring on atom chains of odd numbers by 
C^(Ai, A 2 ). We note that 

7 d< v / 3£#4' ) (Ai,A 2 ). (24) 

The projection onto the X2-axis of the spring in v 7 -direction is an interval J 
of length 75. Counting broken springs, it is elementary to see that the energy 
contribution Jj £ l _ £ ^ )X j W/\{Vy £ ) of the part of these broken triangles that lies 
in the stripe (0, 1) x J is bounded from below by 

2e(l + P( 7 ))^, (25) 

where P( 7 ) is the projection coefficient from (JTj) satisfying P(l) = \ and in 
particular P( 7 ) = if and only if 7 = On the other hand, the energy within 
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stripes (0, 1) x J' when J' is the projection of an arbitrary broken triangle is still 
bounded from below by 2e(3 fl . 

Now let Aj, i = 1, . . . ,M £ , denote all triangles A in C £ ~ R such that 

there exists ^ ^ with (0, 1) x {x^} intersecting with the interior of A. The 
numbering shall be chosen so as to satisfy x% < . . . < x 2 . As 1 — = 0(e), 
there exists a constant C > such that x% + ^ — x£ < Ce, i = 1, . . . , M £ — 1. 

We define the subset {^2 }i=i,-JV" £ of {i^J^l...^ such that = x% for 
a j = l,...N e if and only if — mA i+1 | > 0. According to our previous 
considerations, if 1% is the projection of : = Uj^i w 7 (A^., A ij+1 ) onto the 
^2-axis, then 

|4 7 |< T £#C«. (26) 
As before using (1251) and (126]) we see that the total energy is greater or equal to 

#4?2e(l + P( 7 ))/^ + |P \ PjJ^- + 0(e) 
= + 2#C«£P( 7 )/^ + 2$C$e$* - |4 |— + 0(e) 

ry 7 7 7 ry 

>^ + 2#4 1 7 )eP( 7 )^ + 0(e), 

7 

and so #(3^ = 0(1). As every A e C^; is in at most two different Ci 7 (A i:P Aj. + i), 
this also yields Ef=i #$?(A< ., A ij+1 ) = 0(1). 



Applying (1211) we find that 



- jdi c 

\m&. - m Ai 



0(l) = E#C«(A ij ,A, j+I )>i:^>iE 
j=l j=i v J=1 



for a constant c > 0, when Jj = 7 ~ 1 |?tia ! — ^a !+ J- This concludes the proof. □ 
Remark. We recall that the last claim in Lemma 15.11 does not hold if v 7 is not 
unique (7 = ^) as the crack can take a serrated course. Indeed, if P( 7 ) vanishes, 
we cannot conlude that ftC^ = 0(1) in the above proof. 

The above Lemma 15.11 shows that for a sequence of almost minimizers (y £ ) 
satisfying fl9]), the number j^C £ of largely deformed triangles is bounded indepen- 
dently of e for a < a ciit , while in the supercritical case there are two subsets 

QW : = { x e tt £ : < x x < p £ - ce + (v 7 • ei)x 2 } , 
:= {x G Q £ : p £ + ce + (v 7 ■ ei)x 2 < x± < 1} , 

c > independent of e and p e to be chosen appropriately, such that the number 
of triangles in C £ intersecting Q £ U£l £ 2 ^ is bounded uniformly in e. The following 
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lemma shows that broken triangles in these sets can be 'healed'. In order to treat 
these cases simultaneously in the following we will call these sets the 'good set' 



^good - S (1) (2) 



Q £ for a < a cr it and 

U for a > ocra. 



Lemma 5.2 There exists y £ G W 1 ' 00 (Q goo d, ffi 2 ) with Vy £ bounded in L°°(f2g OOC i) 
uniformly in e such that 

\{x e Q good : y e {x) ^ y £ (x)}\ = 0(e 2 ) 

and 

f dist 2 (Vy e (x), SO (2)) dx<C I dist 2 (Vy e , SO (2)) dx. 

Proof. For notational convenience we drop the subscript e in the following proof. 
We can partition the area covered by the (closed) triangles in C intersecting Q goo d 
into connected components Ci, . . . , Cn such that 

(J A = C 1 U...UCV, 

AeC:AnQ good ^0 

where iV is bounded uniformly in e. Then the maximal diameter of each sets Cj 
is bounded by a term 0(e). For each i, the largest connected component D{ of 
the complement Q goo d \ Ci lying in the same component of Q goo d is unique (with 
area of the order 1 while all the other components of the complement are of size 
0(e 2 )). Let Vi be the union of triangles whose interior is contained in Di that 
touch the boundary of Cj. 

We now proceed to define y by modifying y on all the triangles not contained 
in Di, successively for i — 1, . . . , N. For each i this modification is done iteratively 
on triangles A which share at least one side with a triangle that has been modified 
previously or with a triangle lying in Vi in such a way that y is continuous along 
such sides and y\& is affine and minimizes dist((y)A, SO (2)). 

In order to estimate dist(Vy, SO (2)) we recall the following geometric rigidity 
result proved in [21] : If U C M. d is a (connected) Lipschitz domain, then there 
exists a constant C = C{U) such that for any / e H 1 ^, ~R d ) there is a rotation 
R G SO(d) with 



\Vf(x)- R\ 2 dx <C / dist 2 (V f(x),SO(d))dx. (28) 
u Ju 

The constant C(U) is invariant under rescaling of the domain. For later use we 
mention that if dist 2 (V/(x), SO(d)) is equiintegrable, then R can be chosen in 
such a way that also |V/(x) — R\ 2 is equiintegrable, cf. [25] . 



26 



Consider a single step in the modification process, when y is modified to y 
on A, and let U be the union of triangles that have been modified previously or 
lie in Vi. By the geometric rigidity estimate (128jl . there is a rotation R £ 50(2) 
such that (128j) holds for f — y. Since is piecewise constant, this means 

\(y)^'- R \ 2 < C Y1 dist 2 ((y) A ',50(2)). 

A'cU A'cU 

It is not hard to see that there exists an extension w of y from U to U U A such 
that 

|H A -i?| 2 <0 ^ |(y)A'-i?| 2 . 

A'ct/ 

(If there is only one side of A on the boundary of U, say adjacent to A' C U, 
then one can take w with (w)a = (v)a'- If at least two sides, say in vi and 
v 2 direction, are shared by triangles Ai, A 2 C U, respectively, then these sides 
have a common corner and the unique extension w satisfies (w)a v ^ = {^A^i — 
Rvi + ((y)Ai — R)vi, i — 1, 2.) Now by construction of y on A we see that 



dist 2 ((y) A ,50(2))<0 Y, \(v)a>-R\ 



2 



A'CU 

and so 



dist 2 (Vy(x), 50(2)) dx<C dist 2 (Vy(x), 50(2)) dx. 

UUA JU 

Iterating this estimate we finally arrive at 

/ dist 2 (Vy(x), 50(2)) dx<C f dist 2 (Vy, 50(2)) dx. 

Here the constant C can be chosen independently of e. This is due to the facts 
that the number of modification steps is bounded uniformly in e and - after 
rescaling the shapes U with | - there is also only a uniformly bounded number 
of shapes U involved in the previous rigidity estimates. Moreover, each triangle 
is covered by no more than three of the sets Vi. 

The uniform boundedness of the number of modification steps also shows that 
\{x £ figood : y(x) 7^ y( x )}\ — 0(e 2 ) and, by definition of C and construction of 
y, that ||Vy||^ ( n good) = 0(1). □ 

Note that up to a set of small size y £ satisfies the same boundary conditions 
as y £ on the lateral boundary. More precisely, there are C (0,1), \Ts | = 
0(e), % — 1,2, such that y £ and y £ coincide on (0,e) x ((0, 1) \ r^ 1 ^) and (/ — 
e, I) x ((0, 1) \ T £ ' ). With these boundary conditions and the geometric rigidity 
estimate (|28|) we can now derive strong convergence results for y £ and even the 
corresponding rescaled displacement u £ = -^=(y £ — id) on f2 goo d- We first consider 
the supercritical case. 
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Lemma 5.3 If a > a cl -it, then - possibly after a rotation by ir on the components 
o/r2 gooc i - there exist sequences s e ,t e G M such that 

\\u £ - (0,s e )\\ H1{n w ) + \\u £ - (al,t £ )\\ H1(Q (2) ) -> 0. 

Proof. We again drop the subscript e. Applying the geometric rigidity estimate 
( 1281) to and to fi (2) , we obtain rotations , R {2) G SO (2) such that 

WVy-R^W^^ <C\\dist(Vy,SO(2))\\ L2{ ^ )y i = 1,2. 

Here C can be chosen independently of e as all the possible shapes of fiW are 
related through bi-Lipschitzian homeomorphisms with Lipschitz constants of both 
the homeomorphism itself and its inverse bounded uniformly in e, cf. [21]. Now 
using that Vf/ is uniformly bounded in L°°, we obtain from Lemmas 15.21 and l 3.5( i) 



J2\\Vy-R®\\ 2 L2m <C j dist 2 (Vy,SO(2))dx 

<C / dist 2 (Vy,0(2)) + x (Vy)rfx 

<c/ H/ A , x (Vy)dx. 

JQ^aW L-^ A 



^good\UASC E ^ 

But, as seen before, 

4 f 2B V 
-=- / ^ A , x (Vy) dx < £*(y) - JL\p\ = 0(e) 

Vie Jn good \\j AeCE a 7 

by Lemma 15.11 and so 

Elivs = o(*»). 



i=l 



By Poincare's inequality we then deduce that there are £W G IR 2 such that 



2 



^ 11^-^.-^11^^ = 0(e). (29) 



We extend y as an H 1 function from Q £ to fiW ( as defined in Theorem 12.31) . 
i = 1,2, such that ([29]) still holds andyi(0,x 2 ) = for x 2 G (0, l)\ri 1} , y~i(l,x 2 ) = 

(2) 

1(1 + a £ ) for X2 G (0, l)\Ie . The trace theorem for Sobolev functions with x\ — 
or xi = / according to % — 1 and z = 2, respectively, gives 



2 



=i 
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In particular, setting £W = and = — la e e±, the first components 
satisfy 



£>i - - i2g " -CnL 2((0 , 1)V if>) = (30) 

i=l 

But then also the constant function 

lag = ^ - - izg (■ - - c?) - (x a - - 4? • -cT) 

is of order e in L 2 ((0, |) \rW) and thus < Ce. An elementary argument now 
yields 

\R® - Id| = 0(e) or \R® + Id| = 0(e) 

and, possibly after rotating by tt, we may assume that \R® -Id| = 0(e). 
Returning to (f3"Uj) and (I29p . it now follows that |£i | = 0(e) and then 

(O.CS^IIhi^ 1 )) + 11^- ( a/ 'd 2) ) 11^1^(2)) = 0(y/e). 

□ 

Strong convergence in the subcritical case can be shown along the lines of the 
proofs of the main linearization results in [30] and [31]. We include a simplified 
proof adapted to the present situation here for the sake of completeness. 

Lemma 5.4 If a < a cr it, then there is a sequence s £ G R such that 

\\u £ -(0,s £ )-F a .\\ Hl{ngood) ^0. 



where F a 



a 
-f 



Proof. We again drop subscripts e if no confusion arises. With the help of the 
geometric rigidity estimate ( 1281) we find by arguing as in the proof of Lemma 15.31 
that 

II - R\\% m <C f W A , x (Vy) dx = 0(e) 

for a suitable rotation R e 50(2) with 

\R±Id\ = 0(y/e) (31) 

and 

\\y±id-C\\m<fl e ) = 0(Ve) 
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for some (el 2 with £1 = 0(y/e) and thus, due to the boundary conditions, 

||«- (0^2)11^(0.) = 0(1). 

In particular, u £ — (Ce)2 e 2 converges - up to passing to a subsequence - weakly. It 
now suffices to prove that ||e(u e ) — F a \\ L 2(Q e \ — > 0, where e(u) = ( - Vtl ' > 2 +v " denotes 
the symmetrized gradient, for then the assertion follows from Korn's inequality. 

To this end, we let V £ (F) = ±W A (Id+ y/eF) and V e , x (F) = V £ (F) + ±x(Id + 
y/eF), so that V £>X (F) -»■ ±D 2 V A (Id)[F, F] = \Q{F) uniformly on compact 
subsets of R 2x2 . Then by frame indifference (see Lemma [3.ip 

W A , x (Id + y/eF) = W A , X U (Id + VlF) T (ld + ViF)) 

/ rr \ (32) 

/ f t + F 1 \ 

with /(F) = v/(Id + F) T (Id + F)-Id-^t£, so that \f(F)\ < C min{|F|, |F| 2 }. 
Then by Lemma l3.5( i) and (1321) x satisfies 

V £ , x (^f^ + ^/(^)) > - e dist 2 (Id + ^ 0(2)) + i x (Id + ^) 

> -dist 2 (Id + v^F,SO(2)) 

, 2 

C 



> - 



Id + v / iF) T (Id + y/eF) - Id 
T + F 1 



^F) 



(33) 



In the sequel we set A £ (F) = + ±f(y/eF). Choose convex functions 

ipk : R 2x2 — > R with linear growth at infinity such that ipi < ip2 < ■ ■ ■ and 
ipk(F) — >• IQ(-F') uniformly on compact subsets of M 2x2 . The previous quadratic 
estimate on \4 jX (y4 £ (F)) from below and the fact that V EyX — > \Q uniformly on 
compacts then shows that we can also choose 5 > and a sequence — > oo such 
that 

K, x (A £ (F)) - 5 X{ i MF )i>r k} \MF)\ 2 > ^ fc (A £ (F)) - 1 



whenever £ (depending on fc) is sufficiently small. 
With (1321) we now obtain that 



- / W AtX (y)dx= [ V e>x (A £ (Vu)) dx 

> / i/; k (A e (Vu)) dx + 5 / X{|A £ (v«)|>r fc }|^4 e (Vn)| 2 rfx - 
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As ipk has linear growth at infinity and -^f(y/eVu e ) < Cmin{|Vw E |, a/eI Vw e | 2 }, 
Vu £ bounded in L 2 , by splitting the integration into two parts according to 
|Vu £ | < M or |Vw e | > M and eventually sending M to infinity, we find 

liminf / ^ (A £ (Vu s )) dx = liminf / ipk(e(u e )) dx. 

When u £ — (Ce)2 e 2 ^ u in H 1 , by Theorem 12.11 it then follows that 
ala 2 ,. 4 f Tr , . ._ .. , 

> liminf — / X{dist(x,an)>fc-i}^fc (e(« 6 )) rfx 
v 3 7n 

+ hmsup^= / X{\MVu,)\>r k }\A £ {Vu £ )\ dx-—=-. 
e^o Vo ./n s v ok 

Using that by convexity of ipk the first term on the right hand side is lower 
semicontinuous in Vu £ and that X{dist(-,9n)>fc- 1 }' ? / ; A; ~> \Q monotonically, we finally 
find by letting k — > oo 

V3 v3 ^ (34) 

+ limlimsup— = / X{\A e {Vu e )\>r k }\A E {Vu £ )\ 2 dx. 



A slicing and convexity argument similar to (fT4l) now shows that J n Q(e(u>)) > 
for all w E H 1 subject to u>i(0, x 2 ) = and W\(l, x 2 ) = al and thus 

limlimsup— = / X{\A £ (Vu £ )\>r k }\A £ {Vu £ )\ 2 dx = 0, 
fc^oo e ^ ^3 

or, in other words, |v4 e (Vu e )| 2 is equiintegrable. By the estimate |V^ iX (F)| = 
|iW AlX (Id + y/eF)\ < C(l + |F| 2 ), (|HHD shows that also 

^dist 2 (Vy £ ,SO(2)) < V £ JA £ (Vu £ )) 

is equiintegrable, so that by the discussion following Equation f|28|) in fact we 
may assume that ^||Vf/ e — R\\ 2 L 2(q £ ) is equiintegrable, too, and \R— Id| = 0(y/e) 
by (13TT) . But then also |Vu £ | 2 is equiintegrable and this together with (134j) yields 

lim4= / Q(e(« e )) = 4= / Q(e(u)) = 
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For some 5 > small enough we finally obtain that 
^ ^[Q(F«)dx 



V3 V3 

= inf(^= f Q(e(w)) - 5\e(w) - F a \ 2 dx : 
I v 3 Jn 

w G H 1 (Q),w(0,x 2 ) = 0,w(l,x 2 ) = al 
<liminf^= I Q{e{u £ )) - 6\e(u e ) - F a \ 2 dx 



ala 2 

- - dli 



— - 51imsixp ||e(u e ) - F a \\ 2 L 2 (Ci 



and therefore lim^o ||e(u E ) — F a \\ 2 L2 ^ n ^ = indeed. □ 

After all these preparatory lemmas, the proof of our main limiting result 
Theorem 12.31 is now straightforward. 

Proof of Theorem \2.3[ Choose p e as in ( |27l) and s e , respectively, s e and t e , as in 
Lemmas 15. 4[ respectively, 15.31 By Lemmas 15.41 and 15.31 u E can be extended as an 
if 1 -function from Q £ to Q, respectively, from fie to Q^ l \ % = 1,2, such that still 

\\u £ -(0,s £ )-F a -\\ H1{n] ^0, 

respectively, 

\\ue ~ (0,s £ )\\ H i {nW) + \\u £ - {al, t £ ) || H i (n(2)) -> 0. 

This completes the proof as by Lemma [5.21 we also still have |{i G E : u £ (x) ^ 
u £ (x)}\=0(e). □ 



6 The limiting variational problem 

Convergence of the variational problems 

We first address the question of T-convergence of S £ . The main idea in the proof 
is a separation of the discrete energy into bulk and surface contributions. The 
treatment of the elastic part draws ideas from [31] and [21]. To derive the crack 
energy, one could use a slicing technique, see, e.g., [12]. Although also possible in 
our framework, we follow a different approach here: We carefully construct crack 
shapes of discrete configurations in an explicit way which allows us to directly 
appeal to lower semicontinuity results for SBV functions in order to derive the 
main energy estimates. 

As a preparation we modify the interpolation y on triangles with large defor- 
mation: We fix a threshold explicitly as R = A\/2 and let C £ be the set of those 
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triangles where |(?/)a| > R- We introduce another interpolation y' which leaves y 
unchanged on A G C £ \ C £ and replaces y on A G C £ by a discontinuous function 
with constant derivative satisfying |(?/)a| < R. In fact, by introducing jumps we 
achieve a release of the elastic energy. Note that y' G SBV(Q). 

More precisely, note that on A G C £ we have |(?/)av| > 2 for at least two 
springs v G V. Indeed, for F G M 2x2 with 32 < \F\ 2 = \F T F\ using ([12]) we 
find f |F T F| 2 = §trace(F T F) 2 < E ve v( v > F T Fw) 2 and so < max v6V |Fv| 4 . 
Hence, \Fv\ > 4 for at least one v G V and at least two springs are elongated 
by a factor larger than 2. For m = 2, 3 let C £)m C C £ be the set of triangles 
where |(^)a v l > 2 holds for exactly m springs v G V. Set v 3 = v 2 — Vi. For 

k = 1, 2, 3 pairwise distinct let hi denote the segment between the centers of 
the sides in Vj and direction and define the set V = hj U hj.- 

We now construct y' G SBV 2 (Q £ ). On A G C £ \ C £ we simply set y' = y. On 
A G C £) 2, assuming |(^)a v «| < 2, we choose such that Vy' assumes the constant 
value (y') A on A with (y') A ^i = (v)a v» and \(y') A v\ = 1 for v G V \ {vj}. 
Moreover, we ask that y' = y at the three vertices and on the side orientated 
in Vj direction. This can and will be done in such a way that y' is continuous 
on int(A) \ hi. We note that the definition of (y') A is unique up to a reflection, 
unless (2/)a v i = 0. We may and will assume that 

dist ((y') A , 50(2)) < dist ((y') A , 0(2) \ 50(2)) . (35) 

For A G C £) 3 we set (y') A = Id and y' = y at the three vertices such that ?/ is 
continuous on int(A) \ Vi. Here, the set Vi can be taken arbitrarily at first. 

We define the interpolation u 1 for the rescaled displacement field by u' = 
-^(y f — id). For future reference we define y' v . as 'variants' of y' satisfying that for 
the jump set in some A G C £ $ we always choose Vj. We note that by construction 
also on an edge [p, q] C d A for A G C £ jumps may occur. There, however, the 
jump height \[u £ ]\ can be bounded by 

\[u' £ ](x)\ < e HV^IL < e ■ ce-3 = Cv ^ (36) 

for a constant c > independent of e and x G [p, g]. This holds since the 
interpolations are continuous at the vertices. 

The following lemma shows that we may pass from u £ to u' e without changing 
the limit. 



Lemma 6.1 If u £ — >■ u in the sense of Definition 2.4 and S(u £ ) is uniformly 



bounded, then Xn £ u' £ — > u in L l {VL), Xn e Vi4 — 1 Vm in L 2 {VL) and % 1 {J u i e ) is 
uniformly bounded. 

Proof. We first note that there is some M > such that 

M 

#C £ < — (37) 
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for all e > 0. To see this, we just recall that every triangle A G C £ provides at 
least the energy einf {^(r) : r > 2}. In fact we may assume that C* = C £ in 
Definition El as for A G C* \ C £ we have |(w e )a| < ^=|(z/ £ )a ~ Id] < ^= and so 

||Vw e ||.L2 (fMLWeA) < ||Vu B ||i3 ( n,\u AeC .A) + l|V^|| L2(UAgc|x _ £A) 

It follows that Xf2 e Vw' e is bounded uniformly in L 2 and, in particular, equiinte- 
grable. Finally, the jump lengths 'H 1 {J u ' e ) are readlily seen to be bounded by 
Cej^C £ < C. But then Ambrosio's compactness Theorem for GSBV [21 Theorem 
2.2] shows that indeed Xn £ Vu' £ — ^ Vw in L 2 (Q). □ 

Proof of Theorem{2J% (i) Let e"^a e = a G [0, oo) for all e. Let u G SW 2 ^) and 
consider a sequence w £ C S'i?l /2 (f2 £ ) with y £ = id + £ ^( a e) converging to 

u in S^BV 2 in the sense of Definition 12.41 We split up the energy into bulk and 
crack part neglecting the contribution ^boundary f rom ^he boundary layers: 

£ e (u s )>e W A ((y £ ) A )+eJ2 W A ((y e ) A ) 

A£C e AeC e 

= ~W f W A (ld + ^Vu' E )+eJ2 E 5^(l(i/.)Av|) (38) 



ApC vev, 

£ |(y e ) A v|>2 



=: £ £ elastic (u e ) + £ £ crack ( 



We note that by contruction of the interpolation u' e we may take the integral over 
Q £ . As both parts separate completely in the limit, we discuss them individually. 

Elastic energy. We first concern ourselves with the elastic part of the energy. We 
recall W A (ld + G) = \Q{G)+uj{G) with sup j^, |F| < pj -> as p -> 0. Let 

X e (x) := X[o,e- 1 /4)(|Vu e (x)|). Note that for F G M 2x2 , r > one has Q(rF) = 
r 2 Q(F). We compute 

S?^ c (u £ ) > -1 jf Xe(x) Qq(V<) + ±w (VeVti^))) <tr. 

The second term of the integral can be bounded by 

/]2 oj(\^£Vu' £ \) 



|v^v<| 



2 



Since Vu' £ is bounded in L 2 and Xe-^^^rjr- converges uniformly to as e — > it 
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follows that 



liminf £ £ elastlc (u £ ) > liminf -= / X e{x)-Q{Vu' £ (x)) dx 

4 r i 

> liminf — / -Q(xn e Xe(x)Vu £ (x)) dx. 
£ ~ >u y o Jo. ^ 

By assumption Xn £ Vw £ — Vu weakly in L 2 . As Xe ~ *■ 1 boundedly in measure 
on f2, it follows xn E Xe^u' £ — ^ it weakly in L 2 (f2). By lower semicontinuity (Q 
is convex by Lemma I3.2j) we conclude recalling that Q only depends on the 
symmetric part of the gradient: 



lim inf £, 



elastic , 



u e ) >-^=J^ -Q(e(u(x))) dx. 



Crack energy. By construction the functions u' £ have jumps on destroyed triangles 
A G C £ . We now write the energy of such a triangle in terms of the jump height 
[u] = u + — u~ . We first concern ourselves with some A G C e $. For the variant 
u' e Vi , i = 1,2,3 we consider the springs in Vj,Vfc direction for j, k ^ i. We see 
that the two parts of the jump set, h Vj , h Vk do not overlap in the projection onto 
the one- dimensional hyperplanes spanned by vf and v^- , respectively Thus, we 
compute 

e%)& Vj = e(y' £ ) A v,- + [y' £>Vi \h Vh = ev, + v^K^k* > ( 39 ) 

where [u' £ yJh Vfc denotes the jump height on the set h Vk . Here and in the following 
equations, the same holds true if we interchange the roles of j and k. We claim 
that 



V,- > £ 4 



(40) 



Indeed, for l-^fw^yj^l < e * this is clear since \(y £ ) A Vj\ > 2. Otherwise, 
applying ( l39l) we compute for e small enough: 



> 



it, 



,ViJ/lv 



- 1 



> £3 



> £3 



1 \ _ 1 1 

1 — e-i ) e 4 — 1 > £4 



2 + £~ 



+ 1. 



Let p > sufficiently small. Applying Lemma I3.5( iv) there is an increasing 
subadditive function ip p with ip p (r — 1) — p < W(r) for r > 1. We define ip p = 
ip p — p. The monotonicity of ip p and (140]) yield 



W(\(y £ ) AVj \)>r(\(y £ ) A ^\-l)>r 



(41) 
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Now for A G C £> 3 we may estimate the energy as follows: 

WM&)a) = ^^(I0Qav,|) 
1=1 

* ( e ~ 4 lKvW) +r (^iKvXj)} =■ wa i3 ((&)a) , 
i=i 

where k = 1, 2, 3 are pairwise distinct. With u$ = v u > we can also write 

^a, 3 (OQa) = n-iE / ^ ( e "'IKvJl) (|Vi • ^1 + |v fc • dH 1 . 

The factors in front occur since "H 1 (/i Vj ) = § and, letting Vj be a normal of 
h Vj , one has • Vj| — and |z^- • vjt| = Consequently, defining </>?(r, i>) = 
^ p (r) (|vj ■ i/| + |vjfe • u\) and 0f (r, v) = ip p (r) (|vj ■ i/| + |v* ■ respectively, we 
get 

3 

y/6e J j , nint(A) 

on every A G C ej3 . For A G C £i2 we proceed analogously. Assuming KjQa v «| < 2 
we compute for the springs in Vj, direction (abbreviated by Vj,k) as i n d3H]) 

£(Ve)A Vj,k = e(j/e)A v j,fc + V^K]^ • (42) 

Note that in this case we do not have to take a special variant of u' £ into account. 
Repeating the steps (HP]) and (HT|) we find 

^WK^A^D + Wd^AV.I))^^ (e-i|[uX|) =: Wa.2 (0/.)a) ■ 

Noting that |vj • = jv^ • Vi\ = |vj • = and that every of these terms 
occurs twice in the sum of the right hand side of the following formula, it is not 
hard to see that this energy satsifies the same integral representation formula as 
Wa, 3 : 

3 

^A, 2 (a)A) = / ^(^iKvji,^W- 

V3e . =1 J J u , nint(A) 

e,Vi 

Let a > 0. Then for £ sufficiently small the crack energy can be estimated by 
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where E^ udA (y e ) compensates for the extra contribution provided by jumps lying 
on the boundary of some A £ C £ . We will show that this term vanishes in the 
limit. 

Now by construction the <fr P {r, v), % = 1,2,3, are products of a positive, in- 
creasing and concave in function in r and a norm in v. Moreover, u' e and its 
variants converge to u in L 1 with Vu' £ bounded in L 2 and thus equiintegrable. 
By Ambrosio's lower semicontinuity Theorem [2j Theorem 3.7] we obtain 

hminf^ rack ( Me )>i= f J2^~ 1 \M^n)dn 1 -CMp-]im S upE p EUdA (y E ), 

where we used that sup £ 'H 1 ( J u i ) < CM for a constant C > by (J37J). We recall 
that ip p (r) — > (3 for r — > oo. Passing to the limit cr — >■ this yields 

liminf^ crack ( M£ ) > 4= / 2f3Y\vv u \dU 1 -CMp-limsupE p eUdA (y £ ). (43) 

Taking fl36|) and fl37]) into account we compute 

limsup / (e~*|[i4]| ) | < lim CM sup \ \ip p (r) — p\ : r < e~3 . ce\ \ 
e^O J dA v 1 1 J 

= CM p. 

This proves limsup £ \E p udA (y s )\ < CMp for some C > 0. We let p — >• in (|4"3"|) . 
This finishes the proof of (i). 

(ii) The basic tool for the proof of the r-limsup-inequality is a density result 
for SBV functions due to Cortesani and Toader [17]. We suppose W(f2,IR 2 ) is 
the space of all SBV functions u £ SBV(fl,M 2 ) such that J u is a finite, disjoint 
union of segments and u £ W k '°°(^l \ J U ,M 2 ) for all k. Then W(f2,IR 2 ) is dense 
in SBV 2 {Vt, M 2 ) n R 2 ) in the following way: For every u £ SBV 2 (tt, M 2 ) D 

L°°(n,R 2 ) there exists a sequence w e £ W(f2,IR 2 ) such that Utielloo < Halloo anc ^ 

(i) u e — >• it strongly in L 1 (f2,M 2 ), V« e — >■ V« strongly in L 2 (f2,M 2 ), 

(ii) limsup £ ^ 0(z/ u Jrf'H 1 < (j)(v u )d'H l for every upper semicontinuous 
function : S 1 — > [0, oo) satisfying <f)(u) = (j){—v) for every v £ S* 1 . 

In the following a problem arises as we cannot control the boundary values of such 
an approximating sequence. To overcome this difficulty, similar to the implemen- 
tation of boundary values for the discrete configurations, we consider functions 
attaining the boundary data in small stripes at the boundary. 

Let r\ > 0. We first construct a recovery sequence for u £ W(f2) with 
•) = for < X\ < T) and Ui(xi, ■) = al for / — rj < X\ < I. We de- 
fine u £ {x) = u(x) for x £ C £ fl Q and let y e {x) = id + y/eu e (x). Then y e £ A(a E ) 
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for all e (here in the sense of (J4])). By u £ ,u' £ we again denote the interpolations 
on Q £ . 

J u is the union of disjoint segments. Up to considering the translated lattice 
C £ + (£ e , 0) for appropriate £ e — > and passing to a slightly smaller 77, we may 
assume that J u R C £ = 0. Let 5 > and define = {x E J u , \[u](x)\ > 5}. Let 
V £ and V 5 e be the sets of triangles where J u and J*, respectively, cross at least 
one side of the triangle (i.e., J u (J* respectively) and the side have nonempty 
intersection). Then 

< < + C#( J.) (44) 

for a constant C > independent of m G W(f2, IR 2 ) and e, where #(J U ) denotes 
the (smallest) number of disjoint segments whose union gives J u . From now on 
for the local nature of the arguments we may assume that J u consists of one 
segment only. We show 

V 8 e cC £ C V e . (45) 

for e small enough. Let A G V 5 e . We see that, if J* crosses a spring v at point 
x * , say, then a computation similar as in ( )42|) together with Vu G L°° shows 



\(ys)A v| 



-L[u(x*)} + 0(V 



> 4 + 0(1) (46) 



Thus, A G C £ for £ small enough. On the other hand, if we assume A G' T> e , then 
for at least two springs v G V we have K&Qa v I < 1 1 Id + y/sVuW^ < 2 for e small 
enough leading to A ^ C e . In particular, it is not hard to see that C £ = C £t 2- 
We claim that 

|«) A |=0(1) iorA^V £ \C £ . (47) 
This is clear for A ^ T> £ as Vw G L°°. For A 6 C £ = C e2 there is a v G V such 
that (y' e )A v = (He) a v = v + 0(v / i)- By Lemma l3ToT i) and (J35l) we get a rotation 
i? £ G 50(2) such that 

\R £ ~ (y' £ )A\ 2 = dist 2 ((^) A ,SO(2)) = dist 2 ((^) A ,0(2)) < CW A ((y' £ ) A ) = 0(e). 

This yields \(y' £ ) A - Id| = 0(y/e) and thus |«) A | = 0(1). 

We note that Xn s u £ — > u in L 1 as u and thus every u £ is bounded uniformly 
in L°° and, tt being smooth away from J u , u £ — >■ « uniformly on f2 £ \ [Jagd ^> 
where | [Jag© ^1 — Letting C* = V £ this shows that u £ — > u in the sense 
of Definition [231 recalling flUD and the fact that \{u £ )a\ = 0(1) for A V £ . We 
next establish an even stronger convergence of the derivatives. Consider Vu £ on 
triangles in C £ \T> £ . As 11 is smooth there, the oscillation on such a triangle, 
osc^(Vu) := sup {||Vu(x) — Vm(x')|| 00 ,x,x' G A}, tends uniformly to zero (i.e., 
not depending on the choice of the triangle). We thus obtain 

/ ||V« e - VuH 2 , < / (osc £ A (Vu)) 2 ^0 
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for e — > 0, so that even Xn e \u AeT> aVm £ — > Vw strongly in L 2 (Q). Note that in 
fact xn e Vw £ — ► Vw in L 2 (f2). Indeed, on the set of broken triangles we get 

|V< - Vu| 2 < I |V< - Vu| 2 + I |V< - Vu| 2 

<C#C £ e 2 + C#(D £ \C e )e 2 -(^) 2 
Using flUD, (USD and #£> £ < Ce" 1 this leads to 

limsup / |V< - Vu| 2 < limsup C#(£> £ \ Pf)e < CU X (J U \ 4) 



and letting <5 — ^ yields the claim. 

We now split up the energy in bulk and crack parts 



S £ (u £ ) = £ £ elastic (w £ ) + £ £ crack (u £ ) + 0(e) 



as defined in ((38]) (Note that the contribution £ ££°undary is of order 0(e)). Re- 
peating the steps in the elastic energy estimate in (i), applying xn £ Vu' e — > Vw 
strongly in L 2 (Q) and Q(F) < C|F| 2 for a constant C > we conclude that 

limsup£ £ elastic (w £ ) = 4= / ^Q(e(u(a;))) dr. (48) 



It is elementary to see that J u crosses 

V3e 



n \j u )^-^ + 0(l) (49) 



springs in v-direction for v G V. It is not restrictive to suppose that consists 
of only one segment. Indeed, by the continuity of [u] in J u we see that J s (u) 
is the union of open intervals. Thus, up to a set of arbitrarily small size, J s (u) 
consists of a finite number of open intervals and then it suffices to consider one of 
these intervals. Consequently, (14"9"|) also holds for the subset J & u replacing K 1 (J U ) 
by H l (Ji). Recalling (JU, the crack energy may be estimated by 



limsup£ £ crack (w £ ) 

< limsup (V(J«) sup \\V(r) : r > 5e~^ + 0(1)) 

+ U\J U \ J S U ) m a x{W(r) : r > 1}) A ^ \„ u . v | + 0(e) 

= ^(J u V + m\^ maxVT(r)^ ^X>« 



vev 
v|. 
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We finally pass to the limit 5 — > to obtain 

limsup£ e crack (u e ) < f V^K-vlcW 1 . (50) 
Jj u vgV V3 

Applying f|4*8|) and ( 150]) then shows that u e is a recovery sequence for u G W(f2) 
with Ui(iCi, ■) = for < X\ < r\ and U\{%\, •) = al for / — 77 < X\ < I. Now 
let u G 5'-BV r2 (r2) fl L°°(f2) satisfying the boundary conditions «i(0, •) = and 
Ul {l, •) = la. Define w" G S£V 2 (ft) by 

{(0, w 2 (2^ - xi,x 2 )) for < xi < 2?7, 

"(li^i - fifc, ^) for 277 < Xl < I - 2 V , 

(la, u 2 (2l — Ir] — xi, x 2 )) for I — 2r\ < x\ <l. 

We approximate u v with a sequence u'" G W(f2) such that limsup^^ £(u n '^) < 
£(u n ). Due to the proof of the density result in [17] the approximating sequence 
can be chosen in such a way that u^- 1 equals u[ in (0, 77) x (0, 1) and (I — 77, /) x (0, 1) 
as u\ is constant in (0, 2rj) x (0, 1) and (I — 2rj, I) x (0, 1). Consequently, we may 
assume that Ui' 3 (xi, •) = for < X\ < 77 and u[ 3 [x\, •) = al for I — rj < x\ < I 
for j large enough. 

Let (ei)i be an arbitrary null sequence. According to the above compu- 
tation for j large enough let (u^')i be a recovery sequence for vP'* satisfying 
u^f G *A(a Ei ) for £j <C 77. In particular, we have Xn s .(w)^ — > u 71 ^ in L 1 and 
Xu s . V(u')^' — >■ Vn 7 ''- 7 in L 2 for z — >• 00 (as before, (£t)ef an d denote the dif- 

ferent interpolations of «^'/). Choose j(z) for z G N such that u* := Uef G A(a £ .) 
satisfies xn e . (^)* — * u v in L 1 , x^ e . V(w')* — >■ Vu^ in L 2 for i — >• 00 and in addition 
^{J^m) < Cn\J u v) and #(J^, j(i) ) < g. Observing that (u)* ^ (u')* only on 
a set of triangles T>* with 

#p* < cn ^ J ^ + c#(j B , JW ) < cnHJuv) + -<- (51) 

^2 ^2 

by (144|) we conclude that xn e . («)* — > u v in the sense of Definition 12.41 Since 

limsup£ £i (ti*) < £{vP) 

i— >oo 

u* is a recovery sequence for u v . Repeating the above arguments we can find a 
sequence rji —> with <C 77; such that u** := ul 1 ^^ G ^4(a £i ) and «** — )■ u in 
the sense of Definition 12.41 (in particular, we choose rji such that 'H 1 (J um ,j{i)) < 
CH}(J u in) < CH}(J U ) holds and therefore the argument in floTj) can be applied). 
Moreover, it is not hard to see that £(u n ) — > £(u) for 77 -> as 'H 1 ( J,^ fl (0, 4?]) x 
(0, 1)) + V}(J u v n (I - 4?7, /) x (0, 1)) -> for 77 -> due to the construction of 
u v . Then 

limsup £i{u**) < £{u). 
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This together with the arbitrariness of {£i)i shows the existence of a recovery 
sequence for u G S BV 2 (Q) C\ L°° (Q) . Finally, we may drop the hypothesis u G L°° 
by applying a truncation argument and taking Q{F) < C\F\ 2 into account. □ 



Analysis of the limiting variational problem 



We finally give the proof of Theorem \2.6\ i.e., we determine the minimizers of the 
limiting functional S directly. An analogous result for isotropic energy functionals 
has been obtained in [27] . We thus do not repeat all the steps of the proof provided 
in [27J but rather concentrate on the additional arguments necessary to handle 
anisotropic surface contributions. 

Proof of Theorem \2.6[ Let v 7 be unique and thus P(j) 7^ 0, P(j) as in ([TJ. We 
first establish a lower bound for the energy £. An elementary computation yields 



^ |v • v\ > |v 7 • v\ + \/3|v 7 • v\ = |v 7 • v\ + 



VGV 



1 y/1-7 2 

±— ex • v ± — v 7 • v 

7 7 



> — | ei -z/|+2P( T )|v 7 -v| 



7 

for v G S 1 . In the first step we used that ^ v6V \{ V7 } v = ±v3v 7 . Thus, we get 

4 f 1 f 2/3 4/3 

8(u)>—= / -Q(e(u(x)))dx+ / — |ei • v u \ + — P(7)|v 7 • v u \ dH l . 
V o Jn z Jj u 7 v 3 

Using the slicing method (see, e.g., [31 Section 3.11]) and applying properties of 
the reduced energy proved in Lemma 13.31 we get 

£(«)>^ ^(eJVu(x 1 ,x 2 )e 1 V0) 2 dx 1 + ^-#S X2 {u)^j dx 2 + £\u), 

where #S X2 denotes the number of jumps on a slice (0, /) x {£2} and 



8\u)= j ^P{ 1 )\v,-v\dU\ 



It is easy to see that the energy of a slice is larger or equal to min {ala?/y/3, 2/3/7} 
leading to inf S > min |a/a 2 /\/3, 2/3/7}. Testing with u cl and u cr depending on 
a, we see that this bound is attained and that u cl and u cr , respectively, is a 
minimizer. It remains to prove uniqueness: 

(i) Let a < a crit and u be a minimizer of S. Then u has no jump on a.e. 
slice (0,/) x {£2} an d satisfies a.e. eJVuei = a by the strict convexity of the 
mapping 1 1— > (t V 0) 2 on (0, 00). Thus, if J u 7^ 0, the crack normal must satisfy 
a.e. v u = ±e2- Taking ^ 7 (m) into account, we then may assume J u = up to 
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an H 1 negligible set, i.e., u G H x {££). We find u\(x\,x 2 ) = ax± + f(x 2 ) a.e. for 
a suitable function /, and the boundary condition Ui(0,x 2 ) = yields / = 
a.e. In particular, ef Vu e 2 = a.e. Applying strict convexity of Q on symmetric 
matrices (Lemma 13. 2h we now observe e^Vw e 2 = — | and e^Vu e2 + e^Vw ei = 
a.e. So the derivative has the form 

Vu(x) = _°aj for a.e. cc. 

Since Q is connected, we conclude u(x) = (0, s) + F a x = u el (x) a.e. 

(ii) Let a > a crit and u be a minimizer of £. We again consider the lower 
bound for the energy £ and now obtain that on a.e. slice (0, I) x {£2} a minimizer 
u has one jump and a.e. efVwei = 0. The arguments in (i) show that Vw is 
antisymmetric a.e. Now the linearized rigidity estimate for SBD functions of 
Chambolle, Giacomini and Ponsiglione [16] yields that there is a Caccioppoli 
partition (Ei) of Q such that 

u ( x ) = + k)XEi and J u = (J d*E h 

i i 

where Aj = —Ai e M 2x2 and 6j G M 2 . (See [3] for the definition and basic 
properties of Caccioppoli partitions.) As £ 7 (w) = 0, we also note that v u _L v 7 
a.e. on J u . Following the arguments in [27], in particular using regularity results 
for boundary curves of sets of finite perimeter and exhausting the sets d*Ei with 
Jordan curves, we find that 

J u = {Jd*Ei C (p,0)+Mv 7 

i 

for some p such that (p, 0) +Mv 7 intersects both segments (0, /) x {0} and (0, /) x 
{/}. We thus obtain that (Ei) consists of only two sets: E\ to the left and E 2 to 
the right of (p, 0) + Kv 7 , say. Due to the boundary conditions we conclude that 
A\ = A 2 = and by = (0, s), b 2 = (al, t) for suitable s, t G R. □ 
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